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Introduction 


This report describes the estimation of mean gravity anomalies from 
Geos-3 altimeter data. The work itself extended from April 15, 1974 to December 
31, 1977. The main effort was concentrated in the last eighteen months when 
sufficient data became available for large scale computations to be carried out. 

This specific report describes results obtained to fulfill the given 
NASA contract written in response to a proposal made by the author. A number 
of other reports have been prepared related to the altimeter data analysis. These 
will be referenced in the text of this report. 


The Altimeter Data 


All data used for this study was transmitted to us on magnetic tape by 
the Wallops Flight Center. Originally, we received some data in the BCD mode 
(Leitao, Purdy and Brooks, 1975), but later this data was discarded in favor of 
data distributed to us in the binary mode (Ibid.). The data that we received in 
binary form was read using a slightly modified version of the "Caltor" program 
(Leitao et. al. , 1975, revised 1976), Serelis (1975, revised 1976). For our 
specific purposes we created from the original tapes a BCD type of tape but with 
certain additional information that was to be used for editing purposes. We should 
note here that the BCD format contains information that represents averages in a 
certain time interval. This average represents a "frame" value where the time 
extent of a frame depends on the data rate used in the data transmissions. In the 
low data rate the frame length is 2. 048 seconds while in the high data rate the frame 
length is 3.277 seconds. In the low data rate frame there are 20 observations 
while in the high data rate there are either 32 observations (mode 2) or 320 
observations (mode 3). 

The additional information obtained from the binary tapes was as follows 

1. The Calculated Averaged Automatic Gain Control Voltage (AGC) for 
the frame. This value is given on the Wallops binary tape. 

2. The square root of the variance of the AGC voltage within the frame 
as determined from the individual AGC voltage measurements within the frame. 

3. The Radar Altimeter Average Attitude/Specular Gate voltage which 
was given on the Wallops binary tape. 

4. The Hi/s value which is a measure of sea state. This value is 
computed for a frame and is given on the Wallops tape (see revisions dated 3/3/76 
of the Leitao et al. report). 



We should note here that the altimeter data was taken in a global mode 
and in an intensive mode, the latter being the more accurate. For our studies we , 
used only the intensive mode data. 

For a more detailed explanation of the tapes created by us, see the 
report by Herbrechtsmeier (1977). 


Altimeter Data Editing 

The actual data type that we are interested in was the separation 
between the reference ellipsoid and the sea surface. This quantity, referred to 
as "Sea Surface Height", is computed as follows: 

(1) SSH = h - (A - R+b) 

where 


SSH is the sea surface height given with respect to a specified 
ellipsoid; 

h is the height of the altimeter above the reference ellipsoid as 
determined from the "known" portion of the satellite; 

A is the altimeter measurement; 

R is the refraction correction; 

b is the a priori altimeter bias, if any. 

If the sea surface and the geoid would coincide, SSH would be the geo id height which 
is of specific interest to us. 

Errors in the sea surface height as computed from (1) could occur if the 
altimeter measurement was incorrect, the orbit of the satellite was in error or the 
altimeter bias was not known. To edit out bad data we initially made plots of the 
sea surface heights in comparison with the approximately known geoid undulations 
over the corresponding arc segment. Visual comparisons enabled the identification 
of potentially bad data. However, it became readily apparent that such a visual 
process could not be effectively done because of the large amounts of data being 
received. Consequently, we developed an automatic editing process' using criteria . 
found reasonable in the visual editing process, and from comments made about the 
data by Wallops personnel. 

For a frame measurement to be accepted for further processing, it 
must meet the following criteria; 
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1. 0.25 < cr ALT < 8.00 meters where a ALT is the "Altitude Measurement 
Standard Deviation" given on the binary tape; 

2. The absolute value of the sea surface height should be less than 

170 meters; 


3. The Automatic Gain Control Voltage should lie within the following 
range: -78 < AGC < -62 dbm; 

4. The square root of the variance of the AGCV should be less than 

2 dbm; 


5. The Average Attitude/Specu'lar Gate voltage should be within the 
following range: 35 < AASG < 999 dbm; 

6. The Hi/ 3 value (sea state measure) should be less than 88 meters. 

Criterion one was established to delete too noisy or too smooth data. 
Criterion two was used to delete data that could be off because of very bad orbits. 
Criterion three checked to see if the AGC voltage fell within normal limits while 
criterion four examines the variation of the AGCV within the frame. Criterion five 
was used to detect data taken over land, ice, etc. (Miller, 1977, p. 15). Criterion 
six was used to see if convergence had been reached in the iterative computatioh- 
of the sea surface height. If it had not been, we felt the data was poor and should ;; - 
be rejected. 


Clearly, the above criteria are somewhat arbitrary. However, they 
appear to be able to detect most, but not all, of the bad data. A more extensive 
editing discussion may be found in Herbrechtsmeier (1977). 

All the intensive mode data received by us was processed through the 
data editing program. The original number of frame averages that we started with 
was 469, 800. Of this number 52, 111 were deleted by the automatic editing process. 
We give in Table 1 information related to the number of measurements deleted by 
one or more of the edit criteria previously described. The ith diagonal elements 
in this table indicates the number of measurements that have been rejected because 
the ith criterion and possibly one or more of the other criteria have been met. * 

The off diagonal i j elements in the table indicate the number rejected because the 
i and j criteria have been met together. The numbers listed in the "independent" 
column indicate the measurements neglected because only the ith criterion had 
been met. 
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Table 1 


Data Rejection by Criteria 
(see text for explanation) 


Criteria 

No. 

1 

2 

3 

4 

5 

6 

Inde- 

pendent 

1 

7455 

1214 

1489 

2883 

3141 

2974 

1888 

2 


7079 

1890 

3980 

4828 

2083 

412 

3 



14145 

3835 

8950 

3408 

3614 

4 




20102 

12539 

4578 

• 4483 

5 





30432 

8951 

6059 

6 






18673 

6474 


In some cases, short one or two data point gaps existed in the data. 

In this case, a simple linear interpolation was used to fill in the gaps, m this 
way 9, 704 interpolated sea surface heights were generated. 

After this analysis it was discovered that certain longitudes were in 
error (see letter dated November 8, 1977 from H. Ray Stanley). To eliminate 
this error (it was found some months before the November 8 letter) all data was 
gone through by hand so that the longitude error and other small errors not detected 
by the automatic edit process could be found. In this process an additional 5, 559 
measurements were deleted, leaving a total of 421, 834 measurements distributed 
in 2003 arc segments. The distribution of these arcs is shown in Figure 1. 

As a last edit, the new data set was examined to find adjacent values 
that differed by more than 3 meters. If such unexpected variations were found, 
they were deleted. A total of 2,540 measurements were so deleted, leaving a 
final data set of 1976 arcs containing 419, 294 "frame" observations. 


The Ocean Tide Model 


The sea surface height should be corrected for the ocean tide at the 
time of measurement. The binary tapes supplied by the Wallops Flight Center 
gave, for each frame average, the ocean tide associated with the Hendershott 
(1972) tidal model. The specific implementation of this spherical harmonic model (to 
degree 9) is described by Goodman (1975, p. 4 - 14). 
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Figure 1 

Distribution of Arcs Selected After Editing 


Recently, a more detailed M 2 tidal model was developed by Schwiderski 
(1977). Tests made by Schwiderski indicate that his model gives better fit to real 
world data than die Hendershott (1972) model. This newer model is defined by a 
set of amplitudes and phases for each l°x 1° block on the earth's surface. Given 
a specific location and time, the M a ocean tide can be obtained through fairly simple 
formulas after appropriate interpolations are carried out. 

Tidal values for the Schwiderski model were compared with approximately 
4300 Hendershott values already given on the altimeter data records. In some cases 
the agreement was good between the models, but differences up to 1. 2 meters were 
found. We consequently decided to compute and use the Schwiderski tide value at 
all alt im eter "frame" data points that had survived our edit process. 


Bias Removal and Models, 


The edited values of the sea surface heights now available are still 
subject to errors caused by an altimeter bias and orbit error. 

The altimeter bias has been discussed by Martin and Butler (1977). In 
their study they estimated that the intensive mode bias (measured minus true) was 
-5.3 ± 0. 21 m. Although we could have adopted this as constant for all our data, 
we chose to solve for the bias and orbit error because of the success of our initial 
efforts (Rummel and Rapp, 1977) in this area. 

The orbit error arises from a number of sources including gravity field 
error, etc. The orbits used for the Geos-3 data were classified by the Wallops 
Flight Center as to their accuracy as shown in Table 2. 


Table 2 

Accuracy Classification of Geos-3 Orbits Used in 
the Reduction of the Altimeter Data 


Code 

Source 

Quality 

Orbit Type 

A 

WFC 

10 M + RMS 

1 DAY ARC 

D 

WFC 

3-10 M + RMS 

1-8 REV ARC 

G 

NSWC 

5 M + RMS 

3- DAY ARC 

J 

NSWC 

3 M + RMS 

3 REV ARC 

N 

GSFC 

1-2 M + RMS 

5 DAY ARC 


WFC = Wallops Flight Center 
NSWC = Naval Surface Weapons Center 
GSFC = Goddard Space Flight Center 
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The majority of the orbits computed in the early period of the 
satellites were from WFC and were of code A. From September 1975 through 
May of 1976 the orbits were of code J; from May through early June 1976 the 
orbits were code N with the remaining orbits code G. It was thus clear that 
we had a diverse set of orbits to work with and that the selection of an 
approximate error model for the bias and radial orbit error would be very 
important. 


Another possible source of error in the determination of the sea 
surface heights were the parameters of the reference ellipsoid. (The values 
used as a reference by the Wallops Flight Center were 6378145 meters for the 
equatorial radius and 1/298.255 for the flattening.) For most gravimetric 
geodesy problems, it is desirable to have quantities (such as geo id undulations) 
referred to a mean or general terrestrial ellipsoid (Heiskanen and Moritz, 1967, 
p. 214). With respect to such an ellipsoid, the global average geoid undulation 
would be zero. Since the adopted ellipsoid parameters may not refer to the 
general terrestrial ellipsoid it is desirable to use a procedure that will yield the 
sea surface heights with respect to the general terrestrial ellipsoid even though 
we may not know specifically the parameters (or at least the equatorial radius) 
of the ellipsoid. 

Thus in developing our error model we wish to consider three error 
sources; altimeter bias, radial orbit error and equatorial radius error. 

In order to adjust our data, we decided to fit (through an approximate 
error model) the sea surface heights defined along an arc segment to the geoid 
heights defined by a set of potential coefficients. We would then find the param- 
eters of the error model that would yield the best agreement between the potential 
coefficient model and the adjusted sea surface heights (which had been corrected 
for tidal effects). 

For simplicity, consider this adjustment for a single arc. We write 
for a single point on the arc: 


(2) N b + BA = N fi + v 

where 


Nb is the sea surface height based on the Geos-3 data from the 
Wallops Flight Center; 

Nr is the geoid undulation implied by a set of potential coefficients; 
A represents parameters of the error model; 
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B is the design matrix or coefficients of the error model 
parameters; 

v is the residual after the adjustment 


The adjusted (bias removed) sea surface height would be: 


( 3 ) 


NaOj — Nb + BA 


For our work the geoid undulations were computed from the following (Rapp, 1971) 
Ahax a i 

( 4 ) Nr = ~ \ (-7) \ <C A » cosmX+Sjj.sinmX) (sin#) 

&—Z n=0 

where: 


GM 
r ' 

a 

Cjjn » 

<p, X 

y 


is the geocentric gravitational constant; 

is the geocentric distance to the point at which N r 

is being computed; 

is an equatorial radius; 

are fully normalized potential coefficients with Cs,o 
and C 4>0 referred to an ellipsoid of a specified 
flattening; 

are the fully normalized Legendre functions; 
are the geocentric latitude and longitude; 
is the normal gravity at <p 


A number of different error models and £ MAX values were used to deter 
mine their effect on the adjusted undulations as given by equation (3), and on the 
fit itself as represented by the root mean square residual. The error models 
tested are shown in Table 3. 


• Table 3 


Form of Error Model Tested 


Model 

Form 

1 

Xi + X a cos 0 + X 3 sin 0 

2 

Xi 

3 

Xi + Xa (t-to) 

4 

X x + Xa (t - to) + X 3 (t - to f 
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Model 1 was suggested'by Rummel (1976) for use with arcs of length greater than 
22° 5. In this model 0 is the spherical distance from .the equator crossing. In 
models 3 and 4 t is the time of observation and t 0 is a reference time. 

Six arcs of various lengths were selected for testing. The arcs 
described in Table 4. 


Table 4 


Arc Information for Model Tests 


Number 

No. of obs. 

Length (km) 

1 

533 

7570 

2 

67 

970 

3 

262 

3600 

4 

106 

1700 

5 

603 

14100 

6 

355 

5230 


The first test was to see the effect of Xhax on the adjusted undulations 
or sea surface height. To find this, the first three arcs in Table 4 were separately 
fitted to the GEM 7 coefficients to degree 16 using error model 1 for arcs 1 and 3 
and model 2 for arc 2. The fit was repeated with A max = 4, 8, and 12 with the 
adjusted undulations being compared to those found from the degree 16 solution. 

The root mean square differences are given in Table 5. 


Table 5 


Root Mean Square Adjusted Sea Surface 
Height Differences vs 4 h ax with GEM 7 


Arc 

Degree 16 Adjustment Minus Adj. to 
Specified X MAX 

4«ax = 12 

-&MAX - 8 

Amax = 4 

i 

±0.6 m 

±3.3 m 

±3. 5 m 

2 

±1.6m 

±1.5 m 

±7.7 m 

3 

±1.2 m 

±3.9 m 

±7.9 m 


We conclude that at the ± 1 meter level the adjusted undulations are sensitive to the 
chosen £ MAX . 
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A second test was run to determine the effect of errors in the GEM. 7- * 
potential coefficients on the adjusted sea surface heights. To do this a modified 
GEM 7 coefficient set was generated by adding to each of the coefficients a number 
dependent on the standard deviation of that coefficient as given by Wagner (1976). 
The previous adjustments were repeated and the adjusted undulations compared. 
These results are shown in Table 6. 


Table 6 

Root Mean Square Adj. Sea Surface Height Differences Between - 
Original GEM 7 and Modified GEM 7 Potential Coefficients 




- £ 

MAX 


Arc 

16 

12 

8 

4 

1 

±1. 4 m 

±1.1 m 

±0.3 m 

±0.2 m 

2 

±2. 3 m 

±1.3 m 

±1. 1 m 

±0. 1 m 

3 ' 

±2. 0 m 

±1.1 m 

±1. 1 m 

±0. 2 m 


We conclude that at the ±1 m level the errors in the potential coefficients reflect 
in the adjusted undulation values when only single arcs are being adjusted. 

We next tested the fits obtained by the different error models using 
the GEM 7 coefficients to degree 16. These results are given in Table 7. 


Table 7 

RMS Residual After Adjustment of Test Arcs 
Using Different Error Models and GEM 7 to Degree 16 



Error ; 

Model 

Arc 

1 

2 

3 

4 

1 

± 2.73 m 

± 3. 67 m 

±3.66 m 

±2. 74 m 

2 

0.75 

2.60 

1.32 

0.75 

3 

1.45 

2.75 

2.20 

1.45 

4 

6.97 

8.12 

8.11 

6.97 

5 

3.40 

4.13 

3.38 

3.36 

6 

2. 81 

3.27 

3.05 

2.81 
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We see that assuming a constant bias (model 2) gives the poorest fit to the data. 
Models 1 and 4 give essentially identical results while model 3 gives somewhat 
poorer results than either 1 or 4. One wouldprefer to choose, from these tests, 
either error model 1 or 4. However, later tests conducted using global data 
sets showed that it was difficult to separate in the adjustment the X 2 and X 3 
unknown model parameters. We consequently decided to adopt model 3 for our 
error model when the arc length was greater than 22 £ 5. Arguments by Rummel 
(1977) suggest that for arc lengths below this, with the use of the GEM 7 
reference field, only a single bias term should be solved for (i.e. error model 2). 
Such a procedure was followed in the actual computations to be described in a 
subsequent section. 


Sea Surface Topography Modeling 

The tests conducted so far assume that the sea surface heights 
correspond to the geoid. In fact, they do not because of the existence of sea 
surface topography caused by prevailing winds, currents, etc. in the oceans. 

We consequently felt that we should try to see if correcting our sea surface 
heights for sea surface topography using available information would be of help. 

To do this, the sea surface topography map given by Lisitzin (1974, 
p. 149) was digitized in such a way (Katsambalos, 1977) that an interpolated value 
could be obtained at a specified point. A mean value of 280 dyn. centimeters was 
subtracted to give the actual sea surface topography with respect to the geoid. 

With this data, error model 1 was used to fit the corrected (for sea 
surface topography) altimeter derived undulations. The results of these fits are 
given in Table 8. 


Table 8 

RMS Residuals After Adjustment With Error Model 3, 
With and Without Sea Surface Topography 


Arc 

With Sea Surface 
Topography 

Without Sea Surface 
Topography 

1 

±2.75 m 

±2.73 m 

2 

0.75 

0.97 

3 

1.45 

1.52 

4 

6.97 

6.99 

5 

3.40 

3.43 

6 

2.81 

2.84 
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We can see that no significant improvement in fit is obtained by incorporating sea 
surface topography estimates in our data. However, slight improvements are 
noted and further tests should be made. For this report we will not consider any 
corrections for sea surface topography. 


The Adjustment Model 

We have discussed a technique for removing bias from the sea surface 
heights using a' simple model such as equation (2) with a specific error model being 
chosen from Table 3. In addition, we can incorporate into our adjustment an 
equation to represent the fact that when two arcs cross, the resultant sea surface 
height should be approximately the same. The two heights will not be identical 
because of sea state and observational noise. We now develop the mathematical 
model considering these ideas using a somewhat different notation than Hummel (1976, 1977). 

We start by writing a matrix form of equation (2), considering an 
unspecified number of arcs. Then (2) becomes: 


(5) N b +. B A = Nr + v 

where 

N b » Nr , v and A are column vectors, and 
B is the observation coefficient matrix. 

Equation (5) is of the following form (Mikhail, 1976, p. Ill): 


(6) A v + BA = f 
where, by analogy with (2), we have: 

(7) A = - I 

f = Nr - N b 


The cross-over model for the i and j th arcs can be written as: 


( 8 ) 


N*j(i) - N a dj(j) - C ij 



where Nadj is the adjusted undulation (or sea surface height) and C y is the cross- 
over discrepancy after the adjustment. Using equation (3), (8) can he written as; 


(9) , N b <i) + Bi Ai - N b (j) - Bj A, = Cu 

'' n f * 

Equation (9) can be written as: 

(10) -C u + B, Ai - B j Aj = N b (j) -N b (i) 
which is of the general form of (6) with: 

(11) C ij = v 

A = -I 

f = N b (j)-N b (i) 

Considering all cross-overs between all arcs considered in the adjustment, we can 
write a general cross-over observation equation as : 


(12) -v c + Bo A « f« 

where 


Be - [B lfJ - .Bj] 


Equation (5) can also be written in the following form: 


(13) 


-Za + B a A = f A 


where 


f A = Nr - N b 


We can combine (12) and (13) into a single matrix equation: 
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We let W A be the weight matrix of the altimeter derived undulations, and W c the 
weight matrix of the cross-over observation equations. Then a solution for A of 
{14) is sought such that: 


(15) v A * W a v a + Vc We v c = a minimum 

Under this specification the general solution of; 


(16) 

Av + B A = f 

is 


(17) 

> 

II 

H 

H> 

where 


(18) 

N = B* (AW" 1 A*) 1 B 

(19) 

t = B t (AW -1 A*)" 1 f 

In our case A = 

-Iso that (18) and (19) become; 

(20) 

N = B l WB 

(21) 

t = B*Wf 


Considering (14) and (16) we also have: 



*B a " 


-fft- 


~W A 0 - 

B = 


; f = 


. w = 



-Be. 

* 

-fc- 


-0 Wc- 


so that equation (17) becomes; 



( 22 ) 


A= [BaW a B a + Be WeBef 1 [B*W*{a+ BcW c f c ] 


Knowing the error model parameters from (22) wfe can determine the adjusted 
undulations from: 


(23) N*>j = N b + BA 
while the undulation residuals would be: 

(24) v=N W j-Nr 


The cross-over discrepancy, after the adjustment, would be given by equation (8). 

If there were no cross-over observations, the normal equations in (22) 
would be block diagonal where the size of the block would correspond to the number 
of unknowns for the error model for the specified arc. The block diagonal nature 
of the normal equations reflects the fact that without the cross-over observations, 
each arc can be adjusted separately. 

The cross-over observation equations destroy the block diagonal nature 
of the normal equations. The resulting normal equations in (22) will have the block 
diagonal terms, and also the off diagonal terms reflecting the cross-over information 
between the different arcs. Since, a priori, there will be no particular pattern to 
the normal equations, the solution of them may be time consuming as the number 
of unknowns to be solved for may vary from 1500 to 4000 or more. 

Although the number of unknowns in the normal equations may be large, 
we also expect many zero elements to appear. This suggests that the order of the 
unknowns in the normal equations be arranged in such a way that the column profile 
(for example) of a set of normal equations be minimized so that most zero elements 
are not stored. Then very efficient computations for the unknowns can be made. 

A procedure for minimizing the column profile in normal equations has been 
described by Snay (1976a, b) for possible use in the solution of normal equations 
arising in the adjustment of triangulation networks. The algorithm and Fortran 
program developed by Snay was used in our altimeter adjustment, enabling a 
large number of arcs and cross-overs to be handled in relatively small amounts 
of computer time on our IBM 370/168 operating under MVS. 
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Adopted Geodetic Constants 


For a number of computations to be performed it will be necessary to 
adopt the best available set of geodetic parameters. For these purposes, we use 
the values given in Moritz (1975): 


a = 6378140 m 

GM = 3. 986005 x 10~ 14 m 3 s“ a 

3a = 1082.63 x 10~ 6 

U = 7.292115 x 10" s rad s -1 

The implied flattening is 1/298. 257 and the equatorial gravity is 978031. 8 mgals. 
At the 0. 1 mgal level there is no difference between the gravity formula of these 
constants and that of the Geodetic Reference System 1967. 


The Actual Altimeter Adjustments 

In order to adjust the altimeter arcs shown in Figure 1, we first 
decided to use the GEM 9 potential coefficients (Lerch et al. , 1977) to degree 20. 
This coefficient set is based only on satellite data with no terrestrial gravity data 
in the solution. In addition, it is given to a higher degree than GEM 7, and tests 
reported by Lerch indicate the GEM 9 is an improved solution with respect to 
GEM 7. 


We developed a procedure to develop the cross-over observation 
equations and discrepancies starting from information given in Brown and 
Schroeder (1977) and kindly provided to us on tape by Brown. With information 
in Brown and Schroeder, a procedure was developed to pass through the edited 
tapes to find all cross-overs in a specified geographic region (which could be 
world wide) . This program would then interpolate to the exact cross-over point 
the sea surface height and develop the observation equation coefficients. A more 
detailed explanation of the crossover determination may be found in Herbrechtsmeier 
(1977). 


To implement the discussion in the previous paragraphs, we decided to 
first perform a primary adjustment using a select set of long arcs that had a good 
geographic distribution. The number of such arcs were chosen to be sufficiently 
large such that good coverage with available data would be obtained, and the 
computer time and storage requirements be within our available means. The 
location of the 854 primary arcs are shown in Figure 2. 
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Figure 2. Location of Arcs Used in the Primary Adjustment 
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After adjusting the primary arcs, adjustments were planned for 
regional areas in which sufficient altimeter data was available. In these regional 
adjustments, the primary arcs were held fixed while the non-primary arcs in the 
region received adjustment using the standard procedure using arc and cross-over 
observation equations. The location of the six regional adjustment areas is shown 
in Figure 3. The designations and locations are given in Table 9. 


Table 9 


Regional Adjustment Locations 


Name 

Location 

Alaska 

<p : 10° to 66° 

X: 200° to 260° 

Calibration 

<p: 10° to 50° 

Xr 260° to 310° 

Indian Ocean 

<p:-70° to 30° 

X: 20° to 120° 

Philippines 

(p :-20° to 40° 

Xs 120° to 150° 

New Zealand 

<p ;-70° to 0° 

X: 120° to 214° 

South America 

<p:-70° to 10° 

X: 214° to 20° 


In carrying out the adjustments the standard deviations given for the 
altimeter measured height were used in computing the Wa matrix. These standard 
deviations ranged from 0.25 to 8.0 meters with an average value being on the order 
of 0. 7 meters. To enforce the crossover constraints (or observation equation) the 
standard deviation of this observation equation was chosen as ± 0. 05 m giving a 
weight of 1/400. In the regional adjustments, a weight of 1/800 was given to those 
cross-overs involving a primary arc and a secondary arc. All observation 
equations were regarded as independent. 

In carrying out the global and regional adjustments it was first necessary, 
in some cases, to carry out an initial adjustment. We then examined the cross-over 
discrepancies after the adjustment. In some cases certain arcs gave excessively 
large (greater than 3. 5 meters) cross-over discrepancies. Such arcs were rejected 
from the data set and the adjustment repeated. Some statistical information 
concerning the global adjustment and each of the regional adjustments is given 
in Table 10. 


From Table 10 it is readily apparent that the adjustment process has 
been successful in the reduction of the a priori cross-over discrepancies. The 
discrepancies after the adjustment are of the expected accuracy of the instrument. 
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Figure 3 


Location of the Six Regional Adjustment Areas 


Table 10 


Information Related to Primary and Regional Adjustments 


S 


Region 

Arcs 

Number of; 

Crossover Discrepancies 

Unknowns 

Obs. 

Cross- 

overs 

a priori 

a posteriori 

Global 

854 

1568 

260789 

11406 

±8.27 m 

±0.56 m 

Alaska 

75 

93 

12739 

1575 

9. 63 

0.60 

Calibration 

345 

, 478 

46617 

28920 

7.77 

0.65 

Indian Ocean 

143 

251 

35587 

4912 

8.77 

0.50 

Philippines 

42 

72 

9634 

662 

13.59 

0.69 

New Zealand 

104 

173 

25486 

3193 

8.27' 

0.60 

South America 

185 

281 

34892 

4981 

9.81 

0.57 


We show in Figure 4 a histogram plot of the cross-over discrepancies 
for the global adjustment. Although tests show these discrepancies are not normally 
distributed, the cross-over discrepancy distribution appears reasonable. 

At this point we have a set of sea surface heights (or geoid heights), 
referred to the general terrestrial ellipsoid, and for which we have removed the 
effects of altimeter bias, orbit error, etc. These adjusted values can now be used 
for making geoid maps, the determination of gravity anomalies, or other such 
dependent quantities. 


The Dependence of the Adjusted Sea Surface Heights on the GEM 9 Potential 
Coefficients 


Since we have fitted, by a first degree polynomial, the altimeter derived 
data, there is some concern about the dependence on GEM 9 in the final result. 

From the results previously reported in Tables 5 and 6, we have seen the dependence 
of the adjusted undulations along single arcs on potential coefficient variations. 
However, what may hold along a single track may not hold when the cross-over 
observation equations are introduced. 

To see the effect of potential coefficient changes on the adjusted values 
of the sea surface heights, we repeated the primary adjustment using a modified 
set of GEM 9 potential coefficients to degree 20. This modified set was determined 
using the following equations; 
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where m represents the standard deviation for the GEM 9 coefficient as given by 
Lereh et al. (1977) and d and e are normal deviates drawn from a normal 
distribution with zero mean and a variance of 1. The fee tor 2 was used simply 
to provide a pessimistic error statement. 

The modified potential coefficients were then compared to the original 
coefficient set. The global undulation difference was ±3. 5 meters with the largest 
difference being 11 meters. 

The complete primary adjustment was redone with this modified GEM 9 
potential coefficient set. The adjusted sea surface heights were then compared to 
the corresponding values found from the original solution. The mean difference 
between the two sets of sea surface heights was 0.07 m while the root mean 
square difference was ±0.85 m. 

We thus see that a global difference of ±3. 5 meters in the potential 
coefficient sets has been reduced by a factor of about four in the adjusted sea surfece 
heights. If we accept the original Lereh et al. accuracy estimates as given, the 
root mean square undulation difference would be about ±1. 7 m and we would expect 
this error would propagate into about ±0. 4 m in the adjusted sea surface heights. 
This error is slightly below the average noise level of our observations. 

We thus conclude that errors in the GEM 9 coefficients will effect our 
final adjusted values, but that any error is reduced by a factor of about 4 and in 
our case, this error will be on the order of ±0.4 m. 


Sea Surface Heights or Geoid Undulations as Derived from the Adjusted Data 

The usual geoid maps can be produced from the adjusted data if we 
assume the geoid and sea surface heights coincide. If we do not make such an 
assumption, sea surface height maps can be constructed that should be associated 
with a certain time period. In our case, we choose to produce such maps noting 
that they refer to a time period from April 1975 through approximately November 1976. 

A procedure for constructing these contour maps is described by 
Kearsley (1977). He has used the regional adjustments to construct 2 m contour 
maps in the areas shown in Figure 5 which has been taken from his report. 

Basically, Kearsley has used a least squares prediction procedure, including 
data noise, to predict a geoid undulation and its accuracy on a uniform grid. 
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Figure 5 

Index of Geoid Maps from Geos~3 Altimetry 
(from Kearsley, 1977) 
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These values were then contoured using a standard contouring algorithm. The data 
noise was taken as the altimeter measurement standard deviation. No error com- 
ponent due to the errors in the GEM 9 coefficients nor errors in the adjustment 
process were incorporated. Consequently, the undulation accuracies may be 
somewhat optimistic. 

For purposes of this report, I give two geoid maps taken from the 
Kearsley report. The first is map III which is off the east coast of the United 
States. This geoid map and its corresponding accuracy map are given in Figures 
6 and 7. No altimeter tracks are plotted on these maps because of the great density 
of the available data. This geoid (or sea surface height) was compared by Kearsley 
to a geoid (Nr, r) derived from the GEM 6 potential coefficients and l°x 1° terrestrial 
gravity data by Rapp and Rummel (1975, Table 5), and to the geoid (Nm,c) by Marsh 
and Chang (1976) who used the GEM 8 potential coefficients and l°x 1°, 15* x 15* 
and 5' x 5’ mean gravity anomalies. The results of these comparisons are shown 
in Table 11. 


Table 11 

Comparisons of the Altimeter Derived Geoid with Ground Truth 
Determined from Potential Coefficients and Terrestrial Gravity Data 



Nalt -Nr,r 

Nalt - Nm,c 

Mean Difference 

- 0. 3 m 

- 1.0 m 

RMS Difference 

±1.2 

± 1.5 

[Variance]^ 

± 1.2 

±1.1 

Number of Comparisons 

337 

335 

Maximum Difference 

5.0 

9.2 

Minimum Difference 

- 3.5 

- 4.6 


From this table we can see that there is a greater systematic difference with the 
Marsh, Chang geoid than with the Rapp, Rummel geoid. The RMS differences and 
square root of the variance is consistent with the accuracy of the ground truth geoid, 
being on the order of ±1. 6 m (Rapp and Rummel, 1975, Table 6). A map of the 
undulation differences (both N r ,r and Nm,c) with the altimeter geoid reveals slopes 
that could be indicative of small errors in the potential coefficient sets used in this 
analysis (Marsh, 1977, private communication). 
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The second area in which we give the geoid map (and its accuracy) 
from Kearsley is in the Philippines area (Area I in Figure 5). This information 
is shown in Figures 8 and 9 where the altimeter tracks are represented by the 
dashed lines. From Figure 8 we can see the rapid undulation variation across 
the Philippine Trench, the Mariana Trench and the Izu Trench. The undulation 
variation between the Philippine Trench and the Mariana and Izu Trench is 
smooth. 


Gravity Anomaly Recovery - Theory 

The main proposed goal of our investigation was the determination of 
mean gravity anomalies from the Geos -3 altimetry data, A number of investiga- 
tions were carried out under this contract to determine possible techniques for 
the anomaly recovery. Three such methods are described in a report for this 
contract by Rummel, Sjoberg and Rapp (1977). These methods involved 1) the use 
(by Rummel) of the inverse Stokes equation or the Molodensky equation to compute 
anomalies from geoid undulations; 2) the application of a spectral analysis approach 
(by Sjoberg) to the anomaly recovery and (3) the use of least squares collocation 
(by Rapp). These methods are detailed in the Rummel, Sjoberg and Rapp report 
with each method being tested with altimeter derived undulations that had been 
determined through the adjustment procedure described in the previous sections 
of this report. 

For the purposes of prediction anomaly recovery, we have chosen the 
method of least squares collocation (Moritz, 1975). For the specific equations for 
the anomaly recovery, we have: 


(26) Ag= C gb (Cj* + D) 1 h 

(27) m 2 g = C gs - C gb (C hh + p) -1 Chg 
where 


Ag; the predicted gravity anomaly (either a point or mean anomaly 
depending on the interpretation of the covariances); 

h; a column vector of the given geoid undulations (or sea surface 
heights); 

C g h : the row vector containing the covariance between the anomaly 
being predicted and the given geoid undulation; 

Chh: the square, symmetric matrix containing the covariances between 
the given geoid undulations; 

D: the error-covariance matrix of the given geoid undulations which 

was taken, in our computations, as a diagonal matrix whose 
elements corresponded to the square of the standard deviation 
of the altimeter measurement; 
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Figure 7 


Geoid Undulation Accuracy in Geos-3 Calibration Area 
From Adjusted Altimeter Data 
Contour Interval: 0. 5 meters 
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Figure 8 


Geoid Undulations in the Philippines 
From Adjusted Altimeter Data 
Contour Interval ; 2 meters 
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Figure 9 


Geoid Undulation Accuracy in the Philippines 
From Adjusted Altimeter Data 
Contour Interval: 0. 5 meters 
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C 86 : the expected mean square value (variance) in a global sample, 
of the anomaly being predicted; 

m g : the predicted standard deviation of the predicted anomaly 


We not only intend to predict the mean anomalies in a block, but the mean geoid 
undulations. To do this we write a similar set of equations to (26) and (27): 


(28) 


N = C Nh (Cun + D) 1 h 


(29) 


m n •— C NN — C N& (C ijh + D) C hN 


In this case • 

N: the predicted geoid undulation; 

Cni> : the row vector containing the covariance, between the undulation 
being predicted and the given geoid undulation; 

Cnn: the expected mean square’ value (variance), in a global sample, 
of the undulation being predicted; 
m.N : the predicted standard deviation of the predicted anomaly 


We note that the matrix to be inverted in (26) and (27) is identical to the one to be 
inverted in (28) and (29). 

The interpretation of the predicted quantity as a point or mean value 
depends on die interpretation of or determination of the covariances that are used 
in equations (26) through (28). In our computations we intend to use the point 
covariance function developed in Tscheming and Rapp (1974) on the basis of a 
hypothetical anomaly degree variance model. If we are predicting point undulations 
and point anomalies, the covariances from subroutine COVA (in Tscheming and 
Rapp) can be used directly. If we are determining mean anomalies or mean 
undulations, the appropriate mean covariance functions must be derived. There 
are a number of ways to do this with the technique chosen here to be one of numerical 
integration. Consider two blocks A a and A a whose dimensions are a x b (for A a ) 
and a' x b' for A s . These blocks are, for now, considered to be located on a plane 
in non-identical positions. Now let C(P,Q) be a point covariance function between 
points P (in A i) and Q ( in A 2 ) . Then we can write (following from equation (7-82) in 
Heiskanen and Moritz, 1967): 


(30) 


Ca !a 3 = , ab w a i b ,V P J J* J b C(x,y,x',y’)dxdy dx'dy* 
. (ab)(a'b') J x =o *'y=o‘'x=o J y=o 
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If we now let: 


a = ndy; b=ndx; a' =ndy'; b’=ndx’ 

we can express (30) in the following numerical integration form: 

* a n n a 

(SI) C Al , a = 4JXII C(x lt yi,x'i,y' 1 ) 

i J 1 J 


The accuracy of this computation will improve as n increases, but so will the 
computation time. One wishes to choose n to be as small as possible being 
consistent with the accuracy objectives of the computations . In the computations 
to be described in this report, we choose n = 10 when dealing with 5° anomalies 
and n = 5 when dealing with l°x 1° anomalies. In addition, we interpreted the 
altimeter measurement to be a point value, although the measurement has a 
finite (but small) footprint. 


Gravity Anomaly Recovery - Test Results 

The implementation of the preceding theory required a number of 
decisions on data selection, reference fields, computational efficiencies, etc. that 
had to be made before production computations could be started. We describe in 
this section some of these computations. 

Grid Interval: 


A number of prediction tests were made to select the 
optimum value of n appearing in equation (31). For 5° anomaly predictions, n 
was chosen as 10 and for 1° predictions, n was chosen as 5. Improvement, caused 
by increasing n, in the anomaly predictions (and accuracy) aould be on the order of 0. 3 
mgals, while for undulation estimation, it would be on the order of 0. 1 meters. 

Covariance Tabulation: 

The computation of the needed covariances from COVA 
(Tscheming and Rapp, 1974) can be expensive when many computer subroutine calls 
are made, such as is done in the numerical integration. To get around this, tables 
were constructed of the needed covariances as a function of the spherical distance 
at 0°05 interval. These tables were stored in the program with linear interpolation 
being used to determine the needed covariance. Such a procedure considerably 
reduced the computer time and caused no errors in the final results at the 0. 1 mgal 
and 0. 1 meter levels. 
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Plane Coordinates; 


The covariance functions are usually given as a function of 
the spherical distance 0 separating the points of interest. This spherical. distance 
is usually computed given the latitude and longitude of the two points involved. This 
requires the solution of the spherical law of cosines and the use of an arc cosine 
subroutine. Since we planned on extensive computations with this procedure, we 
examined the idea of getting away from spherical coordinates and use plane 
coordinates instead. The use of plane coordinates might be justified due to the 
fact that we would be working, for each prediction, with data in a limited part of 
the earth. 


The type of plane coordinate system should be simple, 
should have small length distortion, and should have an x coordinate (positive east) 
that was' independent of latitude changes, and the y coordinate (positive north) 
independent of longitude, so that the numerical integration technique could simply 
be applied. Two such sets, of mapping equations were tested. They were the "Carte 
Parallelogrammatique" (Pearson, 1977, p. 250) and a modified Mercator projection 
where the distance is held true at a specified latitude. Not noting any significant 
differences between the results of the two projections, we elected to use the 
following form of the Carte Parallelogrammatique projection; 


(32) 

x = R cos <p o (A - A 0 ) 

(33) 

y = R (<p - <po) 


where E is a mean earth radius (6371 km), <p 0 , A 0 are the coordinates of the midpoint 
of the 5° block being predicted, and <p, A are the coordinates of the point whose x, y 
coordinates are being determined. The spherical distances needed for covariance 
evaluation were found by multiplying the plane coordinate distance by BAD/R.where 
RAD is the radian conversion to degrees factor. 

Limited tests with predictions with plane coordinates gave no 
changes to predicted N values and standard deviation at the 0.02 meter level and at 
less than the 0. 1 mgal level for the anomaly standard deviation. All further com- 
putations used plane coordinates except for a subroutine computing the block variances, 
which was. retained in spherical form. 
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Data Selection: 


Previous experience (Rapp, 1974) and Rummel and Rapp 
(1977) had shown us two facts: 1) only data with about 3° of the edge of a mean 
anomaly block was needed when a high degree reference field on the order of 
12-20 was used; 2) the computational effort in selecting data for a single block 
and the subsequent matrix inversion required a significant amount (about 
2 min. of IBM 370/168 time). Although this time may seem short by itself, 
when one contemplates over 10,000 separate mean anomaly computations, the 
actual time and expense is not reasonable. 

In order to reduce the required computer time, we decided 
to select the data with respect to the 5° equal area block whose mean anomaly was 
being estimated. Specifically, all data within a 3° border surrounding the 5° block 
was selected for possible use. This data was then selected so that all data outside 
a spherical radius of 0° was detected (ip was usually 4° from the center of the 5° 
block). We then said that the maximum matrix that we would invert would be about 
420 x 420. The available data was then selected to be every nth point where n 
was chosen so that the final data set did not exceed 420 (although some tests were 
made up to 475) . We then used this data set to compute the matrix inverse needed 
in equations (27), (28), (29) and (30). Then the appropriate mean covariance 
functions were determined to predict, first the 5° block mean (of Ag and N), and 
then the prediction of all 1° x 1° mean values that fell within the 5° equal area 
block. Thus, if there are s l°x 1° values in the 5° block, we will save s data 
selections and s inversions of a 420 x 420 (roughly) matrix. The time saving is 
obvious. However, we do potentially increase the error in our prediction because 
the data will not be selected optimally for each block being predicted. This is a 
compromise that one must accept for production predictions. 

A number of tests were made comparing the l°x 1° 
anomalies predicted separately, and predicted as part of the l°x 1° anomalies 
within the 5° equal area block. 

Some tests were run that showed the error in the l°x 1° 
anomaly prediction caused by the 5° data selection was on the order of 0. 5 mgals 
reaching 5 mgals in one case (where the 5 mgals was less than the predicted 
standard deviation of the anomaly). Considering the tests made, the error 
attributed to the anomaly prediction caused by the prediction of the l°x 1° anomaly 
with the 5° anomaly data set is probably on the order of 1 to 2 mgals, considerably 
below the expected ± 6 mgal accuracy in the anomalies to be predicted. 

A number of tests were also made specifying the data 
interval along a profile such that every nth point would be used. In average 
coverage areas with respect to the altimeter data, n was 6. In sparse areas n 
could be 2 or 3 (no n below 2 was allowed for stability reasons in the matrix 
inverse). In the densest area, n could reach 19 . In these latter cases, special 
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tests were made to see if a significant improvement in accuracy could be made by 
reducing n (or increasing the data points selected) with the subsequent increase in 
computer time. Comparing the predicted anomalies with terrestrial anomalies on 
an area off the east coast of the United States, the improvement in going from 
an n of 15 to an n of 12 was only 7%, which was considered negligable. These and 
other tests were used to establish the maximum number of data points to be used 
as around 420, provided they were chosen in a manner such that all areas had data, 
if possible. 


Influence of the Reference Field on the Anomaly Prediction: 

The gravity anomalies and geoid undulations predicted from (26) and 
(28) refer to the reference gravity field to which the given geoid undulations (h) and 
the covariances are referred. Usual practice is to give the final values with respect 
to an equipotential ellipsoid model, but a higher degree model such as GEM 9 could 
be used as a reference. Then (26) and (28) would become: 


(32) Ag - C g hR (Chur + D) 1 (h- h R ) + Ag R 

(33) N = C NhR (C toR + Df 1 (h- h R ) + Nr 


where the subscript R indicates a value given with respect to the adopted reference 
field. The needed covariances can be obtained from subroutine COVA previously 
discussed. The reference undulation, h R , can be computed from equation (4). The 
reference gravity anomaly, Ag R , can be computed from: 

•^max jg, A 

(34) AgR = Y (A-l) (~) ^ (Ggn cos m A + S^sinmX) (sin <p) 

r A —3 0=0 


For reference fields to degree 20, (34) is sufficiently accurate for l°x 1° mean values 
when the coordinates of the center point are taken. For the larger 5° equal area block 
the following integrated form of (34) should be used: 


‘'max 


(35) Ag R = 


GM 


r AA(sin<p 2 -sin Qi) 

A 


r*> — P^ s — — 

) (A - 1) [AAC tt , 0 Pn»t> (sin <p) cos (p d <p 
- J J <Pi 


+ y “ [C^b (sin mX 2 -sinmX 1 )- Sg B (cos m X 2 - cos mX^] 


i=i 


J (p 2 _ _ 

P,B (sin <p) cos (p 


Vi 
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A similar equation to (35) follows from equation (4). 

We proposed to test the prediction method using an ellipsoidal model 
(AgR - Nr = 0) and a reference model defined by the GEM 9 coefficients taken to 
degree 20. This test used the adjusted data in the Philippines regional adjustment 
for the prediction of 140 l°x 1° anomalies using the ellipsoidal model and the 
GEM 9 model. The predicted anomalies were compared with a set of terrestrial 
anomalies estimated fromacontour map of Watts (1976). The results are given 
in Table 12. 


Table 12 

l°x 1° Anomaly Comparisons in the 
Philippines Regional Adjustment 



Ellipsoidal Ref. 

GEM 9 Ref. 

Mean Diff. (Alt. - Terr.) 
RMS Diff. (Alt. - Terr. ) 
RMS Alt. Std. Dev. 

2.3 mgals 
±13.4 
± 10.4 

1,5 mgals 
±12.3 
± 7.8 


We see from these values that the influence of the reference field in the final 
predicted values is small with the use of the GEM 9 model yielding slightly better 
results than the ellipsoidal mode. The main difference is the 25% reduction of 
the predicted standard error in going to the GEM 9 field, but the predicted 
standard deviation does not include an error contribution from the potential 
coefficients. Because of these tests we decided to use in all further predictions 
the GEM 9 potential coefficients as the reference field. 


Mean Anomaly and Undulation Production Results 

Using the geo id undulations (or sea surface heights) obtained from the 
regional adjustments, and using the procedures outlined in the previous paragraphs, 
301 5° equal area mean anomalies and undulations and 9995 l°x 1° mean anomalies 
and undulations were computed. The 5° equal area blocks division used is that 
described in Rapp (1977). This scheme, originally suggested by Kaula, was used 
because recent estimates of the anomalies in these blocks were available from 
terrestrial gravity sources, Rapp (Ibid.). 
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The areas in which the mean values were determined are shown in 
Figure 10. These areas were chosen so that they maximized the available altimeter 
data. In addition, emphasis was given to the estimation of anomalies in areas where 
little or no terrestrial data was available. Such areas are apparent from -Figures'll* 
taken from Rapp (1977), which reflected the status of the available gravity>datacas-^V 
of August 1976. ‘ ' - ^ '.V'-.-' 


We report these results in the following way: 

1. hi Appendix A we give a listing of the 5° mean anomalies and 
undulations, and their predicted accuracy. 

2. hi Appendix B we give maps, containing shoreline plots, where all 
l°x 1° anomalies and their accuracy are shown. 

3. hi Appendix C we give, for a single geographic area, maps showing 

the l°x 1° mean geo id undulations and their accuracy. Only this area was plotted 
for space reasons in this report. - 

* 

4. A separate listing of all the l°x 1° anomalies, undulations and thei] 
accuracy which is available from the author. 


Five Degree Mean Anomaly Comparisons 

- • The predicted 5° equal area anomalies from altimeter data have been ' 
compared to the corresponding terrestrial values as given in Rapp (1977). These 
comparisons have been carried out with several subsets of the 5° anomalies'.-.* . 
These subsets are: ■ 

1. The seven anomalies (287, 342, 343, 402, 403, 465, 466) that lie 
off the east coast of the United States in an area of dense altimeter coverage and 
good terrestrial gravity data. 

2. The 153 anomalies in which the terrestrial standard deviation was 
better than ±6 mgals. 

3. All 301 5° anomalies independent of the accuracy of the terrestrial 

estimates. * = " ' 

The results of these comparisons are given in Table 13. 
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Figure 10 

Location of Areas Where 5° Equal Area and l°x 1° Mean Anomalies 
and Undulations Have Been Computed 
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Table 13 


Comparison of Altimeter Derived 5° Equal Area 
Anomalies with Corresponding Terrestrial Values 

(mgals) 


Quantify 

Subset* 

One 

Two 

Three 

MeanDiff. (GEM 9 - Terr.) 

1.3 

1.5 

1.7 

Mean Diff. (Alt. - Terr. ) 

0.2 

1.2 ' 

1.8 

RMS Diff. (GEM 9 - Terr.) 

6.9 

10.0 

12.3 

RMS Diff. (Alt. - Terr.) 

1.7 

6.6 

11.1 

RMS GEM 9 Anomaly 

22.6 

14.8 ■ 

13.8 

RMS Alt. Anomaly 

23.9 

16.3 

14.8 

RMS Terr. Anomaly 

24.3 

17.8 

15.9 

RMS Terr. Std. Error 

2.3 

3.6 

9.7 

RMS Alt. Std. Error 

2.5 

2.6 

2.5 

Number of Blocks 

7 

153 

301 


* see text for description of subsets 


Considering the comparisons with subset one, we see excellent agree- 
ment of the altimeter derived 5° anomalies with the terrestrial data. The smaller 
root mean square difference (±1.7) mgals with the altimeter anomalies than with 
the GEM 9 anomalies (±6.9 mgals) is a strong measure of the possibility of 
improving the GEM 9 gravity field. 

The differences seen for subset two of the 5° equal area anomaly set 
reflect the somewhat poorer terrestrial gravity data in this area. The root mean 
square difference between the altimeter derived anomalies and the terrestrial data 
(of ± 6. 6 mgals) is consistent with the standard deviations of the altimeter and 
terrestrial anomalies. 

The comparisons with the full 301 value data set are shown for general 
purposes only as some of the terrestrial estimates are based on predictions using 
no l°x 1° terrestrial anomalies within the 5° block. In examining the differences 
of all the anomalies predicted, a number of large discrepancies are seen with the 
terrestrial data. Specifically, we found 18 values where the differences were 
greater than 20 mgals (in absolute value). I give in Table 14 the blocks and 
pertinent information having anomaly differences greater than 30 mgals. 
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Table 14 


5° Blocks in Which the Difference Between the 
Altimeter Derived Anomaly and the Terrestrial Anomaly 
is Greater Than 30 mgals (in Absolute Value) 


Block 

AgALT 

AgTERR , 

Difference 

1363. 

1. 3 ± 2. 5 mgals 

67. 1± 7.8 mgals 

-65. 8 mgals 

1209 

17. 1 ±2.5 

-34. 4 ± 6.3 

51.5 

1269 

11.0 ±2.5 

-35.6 ±10.5 

46.6 

1272 

15. 2 ±2. 4 

-24.1 ±12.0 

39.3 

1332 

8.4 ±2.5 

-26.3 ±10.1 

34.7 

1386 

9.5 ±2.5 

-22.4 ±12.8 

31.9 


In most cases these discrepancies are much larger than expected 
considering the given standard deviations. In the case of block 1363, the positive 
terrestrial anomaly estimate is based on some l°x 1° anomalies on the order of 
100 mgals. Such magnitudes are not in agreement with the altimeter derived values. 
To explain the large discrepancies shown in Table 14, a more detailed examination 
of the terrestrial gravity data in these areas is needed. 


One Degree by One Degree Anomaly Comparisons 

We have compared the predicted l°x 1° anomalies with existing 
terrestrial gravity data in two areas: the area (in 8 5° equal area blocks) off the 
east coast of the United States having dense altimeter coverage and good gravity 
data; and the area (in 20 5° equal area blocks) in the Philippines regional adjust- 
ment area where a new gravity anomaly map by Watts (1976) is available. In 
addition, the anomaly variations in the Philippines area is much greater than in 
the U. S. calibration area. The results of these comparisons are given in Table 15. 
The increased discrepancy in the Philippines Area (over the calibration area) is due 
to the much larger anomaly variations that exist in the Philippines Area. In both 
areas -the predicted standard deviations of the altimeter derived anomalies (on the 
order of ±7 mgals) appears to be consistent with the root mean square differences 
between the altimeter and terrestrial anomalies. 

In addition to the above comparisons, we also compared all the l°x 1° 
altimeter derived anomalies to the anomalies (when available) on the August 76 
tape (Rapp, 1977). We found that there were 5119 altimeter anomalies where 
there were no terrestrial estimates available. In addition, we found a number of 
blocks where the discrepancies were quite large between the two anomaly types. 
Specifically, we found 157 blocks where the absolute anomaly difference was 
greater than 50 mgals, and 10 blocks where the difference was greater than 100 
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Table 15 


Comparison of Altimeter Derived l°x 1° Anomalies 
with Corresponding Terrestrial Values 
(mgals) 



Calibration Area 

Philippines Area 

• Mean Diff. (GEM 9 - Terr. ) 

1.8 

-0.4 

Mean Diff. (Alt. - Terr. ) 

1.3 

0.1 

RMS Diff. (GEM 9 - Terr. ) 

±16.3 

±30.8 

RMS Diff. (Alt. - Terr.) 

± 8.0 

±15.8 ' 

RMS Alt. Anomaly 

26.1 

31.0 

RMS Terr. Anomaly 

28.2 

33.4 

RMS Terr. Std. Error 

12.4 

10.7 

RMS Alt. Std. Error 

8.3 

7.3 

Number of Blocks 

245 

520 


Table 16 


Selected l°x 1° Mean Anomalies Across 
the Philippine and Mariana Trench 
(mgals) 



Philippine Trench 

Mariana Trench 

<P° 


Ag<3EH9 

^6 act 

Ag(TERR) 

Q 


AgGEM9 

Ag A LT 

Ag(TERR) 

11 

wm 

35 

59 ±11 

45 ±10 

14 


7 

39 ±10 

27 ±10 

11 


33 

85 ± 9 

70 ±12 

14 

1 

7 . 

47 ± 7 

52 ±10 

11 

125 

31 

140 ± 9 

153 ±20 

14 

144 

6 

102 ± 6 

92 ±20 

11 

126 

29 

-97 ± 6 

-107 ±25 

14 

145 

6 

1 ± 8 

40 ±15 

11 

127 

27 

21 ± 6 

31 ±10 

14 

146 

5 

-135 ± 7 

-148 ±25 

11. 

128 

25 

37 ± 6 

43 ±10 

14 

147 

5 

- 47 ± 7 

- 38 ±15 

11 

129 

23 

31 ± 6 

38 ±10 

14 

148 

5 

18 ± 7 

22 ±12 


mgals. The largest difference found was 169 mgals at a block whose northwest 
comer was -32° (<p), 181° (X). Here the altimeter anomaly was -16 ±6 mgals 
while the terrestrial estimate was 153 ±19 mgals. A complete list of the 
locations of these large discrepancy blocks is available. 

The success of the 1° predictions in reliable ground truth areas suggests 
that bad or questionable terrestrial data may be screened by comparisons with the 
altimeter derived anomalies. 
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To demonstrate the ability of the altimeter derived anomalies to re- 
produce the rapid variations of the anomaly field in certain cases, I give in 
Table 16, l°x 1° mean anomalies across a portion of the Philippine Trench and 
i portion of the Mariana Trench from GEM 9 potential coefficients, from altimeter 
lata and as estimated from the Watts (1976) map. 

From this table we can see the rapid variation of the anomaly that is 
mplied by the altimeter data and is not seen in the low degree anomalies derived 
from a set of potential coefficients given to degree 20. We also see a good agree- 
ment between the altimeter and terrestrial data. 

Preliminary computations have shown that it is possible to produce 
essentially a continuous anomaly profile from the altimeter data that agrees quite 
veil with estimated values based on terrestrial data. With such a procedure, we 
lave seen anomaly changes of 450 mgals in 130 km. 


Summary and Conclusions 

In this report we have described the methods used to determine mean 
anomalies and mean undulations from the Geos-3 altimeter data available to us by 
She end of September, 1977, without having a complete set of precise orbits. The 
sditing of the data was extensive to remove questionable data, although no filtering 
3f the data was carried out. An adjustment process was carried out to remove orbit 
srror and an altimeter bias. This adjustment was carried out first in a primary 
net, and then in regional nets. This adjustment used a very simple first degree 
polynomial in time for the orbit and bias error (for sufficiently long arcs) and 
Introduced cross-over observation equations. After the adjustment, a set of 
adjusted geoid undulations were available that were used to produce geoid undula- 
tion maps (of about ±1 m accuracy) and mean gravity anomalies and undulations in 
301 5° equal area blocks and 9995 l°x 1° blocks. The predicted accuracy of a 5° 
oloek was ± 3 mgals while it was ± 6 mgals for l°x 1° anomalies. These accuracies 
indicate a significant improvement (by a factor of roughly 2) over the accuracies that 
sould be obtained from the GEM 9 potential coefficients themselves. Values were 
predicted only where the altimeter data was of sufficient density°to yield best results. 
Other areas could also be predicted with somewhat poorer accuracies. 5119 l°x 1° 
anomaly estimates took place in areas where no terrestrial data had previously been 
available. Where terrestrial data was available, comparisons with the altimeter 
derived data were made. For the most part, excellent agreement was found between 
the two sets of anomalies indicating the predicted standard deviations are reasonable. 

In the future, more altimeter data is needed to fill in the gaps obvious 
in Figure 1. When such data is available, additional anomaly estimations can be 
made. Because of the homogeneity of the altimeter data, and. because of its 
distribution, it would appear that anomalies derived from satellite altimeter may 
ae more reliable, as a whole, than those derived from ship board data. 


- 42 - 



We finally note that additional work needs to be done to improve the 
orbit error model, and to investigate the errors caused by sea surface topography. 

The gravity anomalies derived from the altimeter data appear to be 
an excellent source of information for improving our knowledge of the structure of 
the earth. 
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Appendix A 


This appendix contains a listing of the 5° equal area free-air 
anomalies and undulations as derived from Geos-3 altimeter data. These values 
refer to an ellipsoid whose flattening is 1/298. 257 and whose equatorial radius 
is theoretically unknown, being actually the equatorial radius of .the general 
terrestrial ellipsoid. 
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SEQ. 

LAT 

LONG 

ANOM 

S.D. 

UHD 

S.D. 

SEQ 

LAT 

LONG 

ANOM 

S.D. 

TJHD 

S.D. 


134 ' 

60 

55 

212 203 

15.0 

2.9 

12.4 

0.2 

135 

60 

55 

222 212 

14.0 

2.7 

8.3 
- 1 . 1 

0.1 

0.1 


176 

55 

50 

213 205 

23.6 

2.5 

7.2 

0. 1 

177 

55 

50 

221 213 

8.8 

2. 5 


3 78 

55 

50 

229 221 

- 2.7 

2.6 

- 9.4 

0 . 1 

223 

50 

45 

213 206 

. 3.0 

2.5 

—3 . 8 

0*1 


224 

50 

45 

220 213 

- 4.2 

2.5 

- 13.2 

0 . 1 

225 

50 

45 

228 220 

- 14.8 

2.5 

— 23 * 2 

0 . 1 


226 

50 

45 

235 228 

- 6.7 

2.7 

- 21.7 

0 . 1 

275 

45 

40 

217 211 

- 3.9 

2.5 

17 . 1 

0*1 


276 

45 

40 

224 217 

- 12.2 

2.5 

- 26.5 

0 . 1 

277 

45 

40 

231 224 

- 12.2 

2.5 

- 29 . 8 . 

0 * 1 


278 

45 

40 

238 231 

- 7.6 

3.0 

- 28.3 

0.4 

286 

45 

40 

292 285 

0.6 

3.9 

- 29.2 

0.8 


287 

45 

40 

299 292 

- 6.6 

2.5 

- 24.7 

0 . 1 

331 

40 

35 

221 215 

- 21. 1 

2*5 

- 29.7 

0 . 1 


332 

40 

35 

227 221 

- 21.9 

2.5 

- 34.8 

0 . 1 

333 

40 

35 

234 227 

- 24 . 1 

2.5 

— 44. 0 

0. 1 


342 

40 

35 

291 284 

- 23.3 

2.7 

- 45.9 

0 . 1 

343 

40 

35 

297 291 

- 23.4 

2. 4 

- 36.8 

0. 1 


391 

35 

30 

224 218 

- 9 . 1 

2.5 

- 29.8 

0 . 1 

392 

35 

30 

230 224 

- 14 . 1 

2.5 

- 36. 4 

0.1 


393 

35 

30 

236 230 

- 15.0 

2.5 

- 40.4 

0 . 1 

394 

85 

30 

242 236 

- 13.9 

2.6 

- 39.5 

0 . 1 


401 

35 

30 

283 277 

- 3.9 

3.6 

- 33.7 

0.6 

402 

35 

30 

209 283 

- 32.5 

2.5 

— 45 . 8 

0. 1 


403 

35 

30 

295 289 

- 23.9 

2.5 

- 44.8 

0 . 1 

437 

30 

25 

129 124 

18.6 

2.5 

23 . 1 

0 . 1 


438 

30 

25 

135 129 

- 13.8 

2.5 

- 2.7 

0 . 1 

439 

30 

25 

141 135 

20.5 

2.5 

40. 4 

0 . 1 . 


454 

30 

25 

225 219 

- 5.4 

2.6 

- 27.9 

0 . 1 

455 

30 

25 

231 225 

- 12.0 

2. 5 

— 36 . 4 

0 . 1 


456 

30 

25 

236 231 

- 16.3 

2.4 

- 35.7 

0 . 1 

457 

30 

25 

242 236 

— 18.3 

2.5 

- 44. 2 

0 . 1 


453 

30 

25 

248 242 

- 13.7 

2.7 

- 39.4 

0 . 1 

462 

30 

25 

270 264 

- 12.0 

2.9 

- 26.7 

0 . 1 


463 

30 

25 

276 270 

- 3.3 

2.7 

- 26 . 1 

0 . 1 

464 

30 

25 

281 276 

5 . B 

2.6 

- 21.8 

0.2 


465 

30 

25 

287 281 

- 17.8 

2.5 

- 39.8 

0 . 1 

466 

30 

25 

293 287 

- 30.2 

2.5 

- 49.2 

0 . 1 


502 

25 

20 

129 124 

3.4 

2.5 

30.0 

0.1 

503 

25 

20 

134 129 

9. 1 

2.5 

37.2 

0 . 1 


504 

25 

20 

140 134 

14.4 

2.5 

51.5 

0 . 1 

505 

25 

20 

145 140 

14.6 

2.5 

43. 1 

0 . 1 


506 

25 

20 

150 145 

0.9 

2.6 

34.9 

0 . 1 

522 

25 

20 

236 231 

- 16.2 

2.5 

- 42.5 

0 . 1 


523 

25 

20 

242 236 

- 26.6 

2.5 

- 55.5 

0.1 

524 

25 

20 247 242 

- 19.3 

2.5 

- 42.9 

0 . 1 

1 

525 

25 

20 

250 247 

- 9.4 

2.6 

- 39.8 

0.1 

528 

25 

20 

269 263 

- 20.7 

2.6 

- 25.8 

0 . 1 

1 

529 

25 

20 

274 269 

, 6.6 

2.6 

- 16.9 

0.2 

530 

25 

20 

279 274 

-5 . 4 

2.4 

- 19.9 

0 . 1 

00 

531 

25 

20 

285 279 

8 . 1 

2.5 

- 31 . 1 

0 . 1 

532 

25 

20 

290 285 

- 32.5 

2.6 

- 42.4 

0 . 1 

1 

569 

20 

15 

125 120 

16.0 

2.9 

36.4 

0.2 

570 

20 

15 

130 125 

13 . 1 

2.5 

42.0 

0 . 1 


571 

20 

15 

136 130 

5.8 

2.5 

53.8 

0 . 1 

572 

20 

15 

141 136 • 

13.9 

2.4 

49 . 4 

0 . 1 


573 

20 

15 

146 141 

31.2 

2.5 

50.7 

0. 1 

574 

20 

15 

151 146 

- 4.5 

2.4 

39.7 

0 . 1 


591 

20 

15 

240 235 

- 20.4 

2.6 

- 44.3 

0. 1 

592 

20 

15 

245 240 

- 20.7 

2. 5 

- 44.0 

0 . 1 


593 

20 

15 

250 245 

- 12.4 

2.6 

- 37.0 

0 . 1 

594 

20 

15 

256 250 

- 4.3 

2.6 

- 31.0 

0 . 1 


598 

20 

15 

277 271 

13.2 

2.9 

- 7.2 

0 . 1 

599 

20 

15 

282 277 

12.2 

2.5 

- 10.2 

0 . 1 


639 

15 

10 

129 123 

29.3 

2.6 

64.2 

0 . 1 

640 

15 

10 

134 129 

9.9 

2.4 

54. 1 

0 . 1 


641 

15 

10 

139 134 

11.8 

2.4 

55.1 

0 . 1 

642 

15 

10 

144 139 

— 4. 4 

2.5 

51.4 

0.1 


643 

15 

10 

149 144 

- 10.8 

2.5 

45.4 

0 . 1 

710 

10 

5 

132 127 

10.0 

2.5 

62.5 

0 . 1 


711 

10 

5 

137 132 

20.0 

2.5 

64.4 

0 . 1 

712 

10 

5 

142 137 

21.9 

2.4 

63.4 

0 . t 


713 

10 

5 

147 142 

19.3 

2.4 

60.3 

0 . 1 

782 

5 

0 

135 130 

37.9 

2.5 

72.9 

0 . 1 


783 

5 

0 

140 135 

24.7 

2.5 

70.7 

0 . 1 

784 

5 

0 

145 140 

19.2 

2.4 

68.5 

0 . 1 


785 

5 

0 

150 145 

12.3 

2.4 

63.4 

0 . 1 

857 

0 

-5 

150 145 

44.7 

2.5 

75.5 

0 . 1 


858 

0 

-5 

155 150 

21.0 

2.5 

65.9 

0.1 

926 

-5 

-10 

137 132 

21.5 

2.5 

63.0 

0 . 1 


930 

-5 

-10 

157 152 

20 . 1 

2.5 

67.6 

0 . 1 

934 

-5 

-10 

177 172 

- 5.2 

2.4 

37 . 1 

0 . 1 


979 

-10 

-15 

46 41 

- 14 . 1 

2.5 

- 25.4 

0.1 

984 

-10 

-15 

72 67 

- 4.5 

2.5 

- 49.0 

0 . 1 


989 

-10 

-15 

98 93 

- 13.7 

2.5 

- 45.3 

0 . 1 

990 

-10 

-15 

103 98 

- 9.2 

2.5 

- 30.3 

0 . 1 


991 

-10 

-15 

108 103 

1.2 

2.5 

- 14.5 

0 . 1 

992 

-10 

-15 

113 108 

- 1.0 

2.5 

- 1.6 

0 . 1 


993 

-10 

-15 

118 113 

- 14,9 

2.4 

9.2 

0 . 1 

1000 

-10 

-15 

154 149 

11.2 

2.5 

67.3 

0 . 1 


1003 

-10 

-15 

170 165 

11.8 

2.5 

58.4 

0 . 1 

1004 

-10 

-15 

175 170 

22.3 

2.5 

55.6 

0 . i 

*1 
ni S 

1005 

-10 

-15 

180 175 

12.2 

2.5 

46.8 

0 . 1 

1006 

-10 

-15 

185 180 

10.5 

2.5 

38.5 

0.1 

1007 

-10 

-15 

190 185 

3.6 

2.5 

28.6 

0 . 1 

1008 

-10 

-15 

195 190 

- 2.2 

2.5 

19 . 5 . 

0 . 1 

1048 

-15 

-20 

42 37 

- 13.7 

3.0 

- 12.6 

0.3 

1058 

-15 

-20 

94 89 

- 22.4 

2.5 

- 51 . 4 ' 

0. 1 

1059 

-15 

-20 

99 94 

- 22.6 

2.5 

- 46.4 

0 . 1 

1060 

-15 

-20 

104 99 

- 26.0 

2.4 

- 39.0 

o.i 

1061 

-15 

-20 

110 104 

- 25.9 

2.5 

- 31.3 

0 . 1 

1062 

-15 

-20 

115 110 

- 17.7 

2.4 

- 12.8 

*■ 0 . 1 

§>£ 

A *-c) 
ci 

1063 

-15 

-20 

120 115 

- 13.4 

2.7 

1.4 

0 . 1 

1069 

-15 

-20 

151 146 

18.9 

2.9 

61.4 

0.2 

1070 

-15 

-20 

157 151 

19.5 

2.6 

73,7 

0 . 1 

1072 

-15 

-20 

167 162 

17.9 

2.5 

60.4 

0.1 



SEG 

LAT 

LONG 

ANOM 

S.D. 

UNI) 

S.D. 

SEG LAT 

LONG 

ANOM 

S.D. 

UWD 

S.D. 


1073 

-15 -20 

172 

167 

31.4 

2.5 

63.7 

0 . 1 

1074 -15 -20 

177 

172 

32.4 

2.5 

61.2 

0.1 


1075 
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Appendix B 


This appendix contains geographic area maps containing the l°x 1° 
mean free-air anomalies, and (on an adjacent map) the corresponding predicted 
accuracies, as estimated from the Geos-3 altimeter data. These maps have been 
prepared for 30 areas, approximately 30° in latitude extent and 20° in longitude 
extent. An index map for the individual maps is shown in Figure B o • The 
individual maps follow, labeled Figures B x through B 30 with the anomaly map 
given first, followed by the standard deviation map. 

The values given refer to an ellipsoid whose flattening is 1/298.257 
and whose equatorial radius is theoretically unknown, being actually the equatorial 
radius of the general terrestrial ellipsoid. 
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Figure B 0 

Map Showing Areas Where l°x 1° Anomalies and Their Accuracy are Given 
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Appendix C 


This appendix contains a sample plot of the l°x 1° mean undulations 
and their accuracy for one geographic region (Block 3 shown in the index map 
given in Appendix B). 

- v. 

The undulations given refer to an ellipsoid whose flattening is 
1/298. 257 and whose equatorial radius of the given terrestrial ellipsoid. 
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Approximately 2000 GEOS-3 altimeter arcs, supplied by the NASA Wallops Flight Center and 
received by September 1.977, have been analyzed to improve our knowledge of the geoid arid' 
gravity field. The first step in. this process was a complex editing procedure used to 
eliminate bad and unacceptable data. After .this, an adjustment- procedure was used to fit 
the sea surface heights (geoid undulations) implied by the altimeter data to the geoid 
undulations implied by the GEM 9 potential coefficients, in. an adjustment process that 
incorporated cross-over constraints.' The error model used for the fitTwas a one or two 
parameter model which was designed to remove altimeter bias and orbit' error. The adjustment 
process was carried out- in two stages: first in a primary net of global data in 854 arcs; 

second in six regional areas. Typical a priori cross-over discrepancies were +8 meters while 
a posteriori discrepancies were on the order of +0.6 meters. 


The undulations on the adjusted arcs were used to produce geoid maps in 20 regions where the 
data was suitably dense. In addition, the adjusted data was used to .derive 301 5° equal 
area anomalies and 9995 l°xl°- anomalies in areas where the altimeter data was most dense, 
using least squares collocation techniques. The average predicted accuracy of the 5° anomalies 
was +3 mgals, and +6 mgals for the l°xl° values. Comparisons with terrestrial data indicate 
these anomaly accuracies were reasonable. We also emphasized the ability of the altimeter 
data 'to imply rapid anomaly changes of up to 240 mgals in adjacent l°xl° blocks. 


17. Key Words {Suggested by Author(s)) 18. Distribution Statement 


geoid 
gravity anomalies 


GE0S-3 
altimetry 



STAR Category - 

42,46,48 


19 Security Classif. (of this report) 

20. Security Classif. (of this page) i 



Unclassified 

Unclassified 



■ 


Unclassified - unlimited 


For sale by the. National Technical Information Service, Springfield, Virginia 22151 






















